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ABSTRACT 

We present a detailed abundance study based on spectroscopic data obtained with HARPS of two solar-analogue main targets for 
the asteroseismology programme of the CoiJoJ satellite: HD 42618 and HD 43587. The atmospheric parameters and chemical com- 
position are accurately determined through a fully differential analysis with respect to the Sun observed with the same instrumental 
set-up. Several sources of systematic errors largely cancel out with this approach, which allows us to narrow down the 1-cr error bars 
to typically 20 K in elfective temperature, 0.04 dex in surface gravity, and less than 0.05 dex in the elemental abundances. Although 
HD 42618 fulfils many requirements for being classified as a solar twin, its slight deficiency in metals and its possibly younger age 
indicate that, strictly speaking, it does not belong to this class of objects. On the other hand, HD 43587 is slightly more massive and 
evolved. In addition, marked dilferences are found in the amount of lithium present in the photospheres of these two stars, which might 
reveal different mixing properties in their interiors. These results will put tight constraints on the forthcoming theoretical modelling 
of their solar-like oscillations and contribute to increase our knowledge of the fundamental parameters and internal structure of stars 
similar to our Sun. 
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1. Introduction 

Space-borne observatories dedicated to asteroseismology and 
the detection of planetary transits, such as the Kepler (Borucki et 
al. l2010l l or CoRoT (Baglin et al. 20091) satellites, are providing a 
wealth of very high quality photometric observations. Modelling 
these data places very strong constraints on stellar models and is 
currently leading to great advances in our understanding of a va- 
riety of physical processes that take place in stars with widely 
different properties and evolutionary status. In particular, ob- 
servations of solar-like oscillations in late-type, main-sequence 
stars can be used to determine some of their fundamental stel- 
lar parameters with unprecedented accuracy or to gain important 
insights into their internal structure (e.g., Christensen-Dalsgaard 
& Houdek 120101 and references therein). Despite these very 
promising prospects and the dramatic achievements that have al- 
ready been made, a good knowledge of the quantities that are 
not accessible from an asteroseismic analysis (most notably the 
effective temperature and the chemical composition) remains an 
essential ingredient of any successful theoretical modelling and 
is needed for trustworthy inferences about the star properties 
(e.g., Carrier et al. 2005 Basu et al. l20T0^ . In the case of CoRoT, 
in particular, the photometric observations have been supported 
by ground-based spectroscopic ones since the launch date in or- 
der to provide these quantities (see Poretti et al. l2013t . 
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This paper is an effort to obtain such accurate constraints 
and restrict the range of possible seismic models for two 
bright, solar-analogue main targets for the asteroseismology pro- 
gramme of the CoRoT space mission: HD 42618 (HIP 29432; G3 
V, y = 6.84 mag) and HD 43587A (HR 2251, HIP 29860, GJ 
23 1.1 A; GO V, y = 5.70 mag). The latter (hereafter HD 43587 
for simplicity) is the primary of a quadruple system with a pair of 
distant, common-proper-motion components and a closeby M- 
type companion (e.g., Fuhrmann 2011). These two nearby stars 
are located about 20 pc from the Sun and lie towards the anti- 
centre direction of the Galaxy (/ ~ 203° and b — 5°). They have 
been observed by CoRoT during long runs (LRa04/LRa05 for 
HD 42618 and LRa03 for HD 43587) with a typical duration of 
about 150 days. 

The extensive observations of these two stars by CoRoT oSei 
a rare opportunity to examine in some detail the internal prop- 
erties of early, main-sequence G stars and to evaluate to what 
extent they differ from those prevailing in the Sun. This is of 
particular interest in the case of HD 42618, which has been rec- 
ognized as a close solar analogue (e.g., Soubiran & Triaud l2004l l 
and is, as we argue below, the star most similar to our Sun in 
terms of fundamental parameters ever observed by the satellite 
in the course of its asteroseismology programme. On the other 
hand, the binarity of HD 43587 provides supplementary infor- 
mation, such as the dynamical mass, which can be compared 
with the asteroseismic estimate. The full analysis of the CoRoT 
data will be presented for HD 42618 by Barban et al. (in prepa- 
ration) and for HD 43587 by Boumier et al. (in preparation). A 
preliminary report of the results for HD 42618 can be found in 
Barban et al. (i20T3l l. 
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Fig. 1. Differences on a line-to-line basis for each line list between the abundances derived from the Fe i and Fe n lines in the CoRoT 
targets and in the Sun as a function of the lower excitation potential. The results for HD 42618 and HD 43587 are shown in the 
left-hand and right-hand panels, respectively. The mean iron abundance obtained for each line list is shown as a dashed line. The 
line-to-line scatters (in dex) are indicated. 



This paper is organised as follows: Section |2] discusses the 
observations and data reduction. The methods of analysis used 
are described in Sect|3] The results are presented in Sectj4] and 
their reliability is examined in SectjS] Section |6] is devoted to a 
discussion of the abundance pattern of our two targets with re- 
spect to that of the Sun. Finally, our main results and conclusions 
are summarized in Sect|7] 



2. Observations and data reduction 

Three high-resolution, optical spectra of each Co/JoT target were 
acquired during the period December 2010-January 20 11 in the 
framework of the ESO large programme LP185.D-0056 (Poretti 
et al. l2013l l with the HARPS specti'ograph attached to the 3.6-m 
telescope at La Silla Observatory (Chile). In the high-efficiency 
(EGGS) mode, the spectral range covered is 3780-6910 A and 
the resolving power R ~ 2,0 000. The signal-to-noise ratio (S/N) 
at 5815 A, as measured on each exposure, lies in the range 263- 
275 for HD 42618 and 287-392 for HD 43587. The data reduc- 
tion (i.e., bias subtraction, flat-field correction, removal of scat- 
tered light, order extraction, merging of the orders, and wave- 
length calibration) was carried out using dedicated tools devel- 
oped at Brera observatory. The spectra were subsequently put in 
the laboratory rest frame and were continuum-normalised by fit- 
ting low-order cubic spline or Legendre polynomials to the line- 
free regions using standard tasks implemented in the 1RAF0 soft- 
ware. The individual exposures were co-added (weighted by the 
S/N) to create an averaged spectrum, which was subsequently 
used for the abundance analysis. 

A series of five consecutive solar exposures (with a typical 
S/N of ~350) were obtained on 21 December 201 1 with the same 



' IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



instrumental configuration. Owing to the faintness of the observ- 
able asteroids at the time of observation, daytime-sky spectra 
were obtained. The telescope was pointed as close as possible 
to the Sun position to minimize the effect on the line profiles 
of scattering in the Earth's atmosphere (Gray et al. 2000 i. The 
impact on our results of this possible mismatch between our 
skylight spectrum and the true solar one is examined in SectjJ] 
Exactly the same reduction steps as above were applied and the 
weighted average spectrum was used for the solar analysis. 

3. Metliods of analysis 

The atmospheric parameters {T^s, logg, and microturbulence 
^) and abundances of 22 metals were self-consistently deter- 
mined from the spectra using a classical curve-of-growth anal- 
ysis. The Li abundance was determined from spectral synthe- 
sis. In each case, Kurucz plane-parallel atmospheric models and 
the 2010 version of the line analysis software MOOG origi- 
nally developed by Sneden (I1973I I were used. Molecular equi- 
librium was achieved taking into account the most common 
22 molecules. All models were computed with opacity distri- 
bution functions (ODFs) incorporating the solar abundances of 
Grevesse & Sauval (I1998I I and a more comprehensive treat- 
ment of molecules compared to the ODFs computed by Kurucz 
(1 19901 1. More details are provided by Castelli & Kurucz (2004)0 
A length of the convective cell over the pressure scale height, a 
- 1.25, and no overshooting were assumed. Appropriate ODFs 
and Rosseland opacity tables (rounded to the nearest 0.1 dex) 
were used according to the mean Fe abundance derived from the 
spectral analysis (because the mean abundance of the a elements 
is solar within the errors, no models with enhancements of these 
elements were considered). 

Thanks to the similarity in terms of fundamental parameters 
between our targets (especially HD 42618) and the Sun, enforc- 

- See also |http : //wwwuser . oat . ts . astro . it/ castelli/ 1 
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ing a strictly differential analysis is expected to minimize the 
systematic errors that arise either from the data treatment (e.g., 
continuum placement) or from physical effects (e.g., granula- 
tion). This approach is particularly advantageous in that it seri- 
ously mitigates the effects of departures from local thermody- 
namic equilibrium (LTE), which can adversely affect the accu- 
racy of the results. Such non-LTE effects are known to become 
increasingly important with decreasing metallicity. As we show 
below, HD 42618 is slightly more metal poor than the Sun at 
the ~0. 1 dex level. However, the differential effect of departures 
from LTE on the derived iron abundance and parameters appears 
to be negligible at these levels of metal depletion (Lind et al. 
12012 ). To first order, our results relative to the Sun should hence 
not be impacted by the assumption of LTE (the same conclusion 
applies to the Li abundances because the non-LTE corrections 
are expected to be equally very small in the Sun and in our two 
solar analogues; Lind et al. 2009). 

We therefore base our results on the differences on a line-to- 
line basis between the abundances derived in the Co/JoT targets 
and in the Sun (see, e.g., Melendez et al. l2012l for an application 
of this method). In particular, the model parameters (T^ff, \ogg, 
^, and [Fe/H]) were iteratively modified until all the following 
conditions were simultaneously fulfilled: (1) the differences in 
the Fe i abundances exhibit no trend with lower excitation poten- 
tial or reduced equivalent width (the logarithm of the equivalent 
width divided by the wavelength of the transition); (2) the mean 
abundance differences found using the Fei and Feii lines are 
identical; and (3) the iron abundance obtained is consistent with 
the model value. The abundance differences were inspected and 
any outlying iron spectral features removed from the analysis 
(this procedure was also used when determining the abundances 
of the other chemical elements). To determine the zero-point 
solar abundances, we adopted Teff - 5777 K, log g = 4.44 dex 
(Cox l2000l l and adjusted the microturbulence until there was no 
dependence between the reduced equivalent width and the Fe i- 
based abundances. We obtain for the Sun a mean value, ^ ~ 0.86 
km s"', as reported in some previous works (e.g., Bensby et al. 
I2003I I. To assess the uncertainties associated to the exact choice 
of the diagnostic lines used and the line atomic data (oscillator 
strengths, broadening parameters), the analysis was carried out 
using four different line lists taken from the literature: Bensby et 
al. (12003]) . C hen et al. (2000), Feltzing & Gonzalez (l200Tll, and 
Reddy et al. (JSOOS). The study of Feltzing & Gonzalez (l200TT l 
concentrated on metal-rich stars and, as a result, the lines used 
are generally weaker than those in the other line lists. Atomic 
lines significantly affected by telluric features were discarded 
from the analysis (the telluric atlas of Hinkle et al. 120001 was 
used). The oxygen abundance could unfortunately not be deter- 
mined because [O i] /16300 is affected by a strong telluric feature 
in our reference solar spectrum. The Fe i and Fe ii abundance dif- 
ferences are shown as a function of the excitation potential (used 
to constrain Teff) in FiglT] The typical line-to-line scatter is about 
0.02 dex. 

The determination of the Li abundance relied on a spectral 
synthesis of Li i /16708. The lithium abundance was determined 
using the accurate laboratory atomic data of Smith et al. (11998I) . 
A solar isotopic ratio ^Li/''Li = 12.18 was adopted (Asplund et 
al.|2009). The van der Waals damping parameters for the lithium 
components were taken from Barklem et al. (12000) . Spectral fea- 
tures of the CN, TiO, MgH, and SiH molecules were considered 
(more details can be found in an upcoming abundance study of 
the red giants in the CoRoT seismology fields; Morel et al., in 
preparation). A knowledge of the total amount of line broaden- 
ing is needed to perform the synthesis. An attempt was therefore 




Fig. 2. Variation for HD 426 1 8 of the fit quality (colour-coded 
as a function of xl) for different combinations of v sin / and ^rt 
in the case of Fe i /16705. 1 . The best fit is found for v sin / = 0.9 
and = 3.7 km s"' (and [Fe/H] = -0.13 dex). 



made to separate the contribution of stellar rotation and radial- 
tangential macroturbulence, ^rt, through fitting three relatively 
unblended Fe i lines, two of which lie in the vicinity of the Li 
doublet (Fe i /16265. 1, 6703.6, and 6705. 1 A). The other free pa- 
rameter was [Fe/H] (for all stars, the iron abundances found are 
identical, within the errors, with the mean values derived from 
the curve-of-growth analysis). Instrumental broadening was de- 
termined from calibration lamps and assumed to be Gaussian 
with a full width at half-maximum (FWHM) of 0.079 A at the 
wavelength of the Li feature. Values in the range 0.4-0.9 and 
3.1^.2 km s"' were inferred for vsin/ and ^rt, respectively. 
Although such macroturbulent velocities are typical of solar-like 
dwarfs (see fig. 3 of Valenti & Fischer 2005 ), there is a strong de- 
generacy in the determination of these two quantities, as shown 
in Fig[2l Furthermore, the same analysis carried out on the so- 
lar spectrum provides a formal v sin / value that is significantly 
lower than the true one (we obtain vsin; = 0.6 and ^rt = 3.8 
km s '). Although this known value is within the range of ac- 
ceptable solutions defined by the strip in the (v sin /-^rt) plane, 
this shows that these two quantities are poorly constrained be- 
cause of degeneracy problems (see also, e.g., Bruntt et al. 120101 ). 
For both stars, only an upper limit on v sin / of about 3 km s"' 
could be established. The Fourier analysis of the spectral lines 
combined with the availability of dedicated libraries of synthetic 
spectra may allow one to separate the contribution of rotation 
and macroturbulence (e.g.. Gray 1984^. Here, exploratory calcu- 
lations of the mean profile Fourier transform suggest, based on 
the relative height and shape of the sidelobes (in particular the 
squeezing of the first one; Gray 2005 ), that macroturbulence is 
the dominant broadening mechanism. This qualitatively supports 
the results of the line-profile fitting described above. Valenti & 
Fischer (2005) obtained vsin/ = 2.1 km s"' for HD 42618 and 
2.0 km s ' for HD 43587 adopting a macroturbulence derived 
from their calibration as a function of the effective temperature. 
As they discussed, however, their values may be overestimated 
for Teff > 5800 K. The difficulties encountered when trying to 
constrain this parameter from the spectroscopic material rou- 
tinely secured for abundance studies is illustrated by the fact 
that Fuhrmann (i2004j estimated vsin/ =1.0 km s"' only for 
HD 426 1 8 by also adopting a calibrated macroturbulent velocity 
(^RT = 3.3 km s"' in that particular case). A low-rotation rate is 
consistent with their low-activity levels (Hall et al. .2009. Gray et 
al. 120031 Isaacson & Fischer: 20T0l Lubin et al. l20T0l ). HD 43587 
shows no evidence for activity cycles (Baliunas et al. 11995) and 
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has even been proposed to go through a Maunder-like minimum 
(Schroder et al. I2012I I. Our synthesis of Li i /16708 is based on 
the best-fitting values, but adopting any combination of v sin / 
and within the statistically acceptable solutions leads to a 
fit quality that is not noticeably different. The iron and lithium 
abundances were adjusted until a satisfactory fit of the blend pri- 
marily formed by Fe i /16707.4 and the Li doublet was achieved. 
A close agreement was found in all cases between the abundance 
yielded by this weak iron line and the mean values found with 
the equivalent widths (EWs). The linear limb-darkening coeffi- 
cients in the R band were taken from Claret (2000l. The fits to 
the lithium feature for the CoRoT stars are shown in Figj3] The 
absolute solar reference lithium abundance we obtain is log e(Li) 
= +0.92 dex. 




6707.4 6707.6 6707.8 6708 
Wavelength (A) 



6708.2 



Fig. 3. Comparison between the observed and synthetic Li i 
/16708 line profiles for HD 42618 (upper panel) and HD 43587 
(bottom panel). In each panel, the observed spectrum is shown 
(filled circles) along with the solar one for comparison (open 
circles). The best-fitting profiles and those assuming no lithium 
are shown as solid and dotted lines, respectively. The parame- 
ters have been derived using the line list of Reddy et al. (I2003I I. 
The vsin/ and macroturbulence values adopted (in km s ') are 
indicated. The vertical tickmarks in the upper panel indicate the 
location of the ^Li components. 



are necessarily accompanied by variations in the other two. We 
therefore also took into account the fact that, for instance, such 
an error in Teff would lead to different values for log g, ^, and ul- 
timately for the abundances. This covariance term was estimated 
by varying T^^fi by its uncertainty, keeping it fixed, and then rede- 
termining log g and ^ as described above. Once again, the analy- 
sis was repeated using this new set of parameters to estimate the 
effect on the abundances. Finally, the line-to-line scatter, crii„e, 
was also taken into account to compute the abundance errors. 

On the other hand, our results are also affected by three main 
sources of systematic errors. As already mentioned in Sectj2] 
our daytime-sky spectrum is not expected to be a perfect match 
to the solar spectrum. First, a sky spectrum is obviously not an 
observation of a point-like source. A potentially more serious 
concern is also that the line profiles may be slightly filled in be- 
cause of scattering in the atmosphere, with changes in the EWs 
amounting to a few per cent (e.g.. Gray et al. 2000 Molaro et 
al. 2008) . To evaluate the impact on our results, we repeated the 
analysis after artificially increasing our EWs by 3%, which can 
be regarded as a generous estimate (e.g., Ramirez et al. 1201 11 1. 
In addition, we carried out the analysis using MARCS models 
(Gustafsson et al. 120081 ). Note that in contrast to ATLAS9 model 
atmospheres, MARCS models use solar abundances based on 
3D hydrodynamical simulations (Grevesse et al. 2007 1. Very mi- 
nor differences are found with respect to the results from Kumcz 
models, which is indeed expected given the close similarity of 
the atmospheric structure for both sets of models in this parame- 
ter range (Gustafsson et al. 2008 1. Finally, estimating the best fit 
to the lithium line is subjective and is especially uncertain in the 
Sun and in HD 42618, which have weak features. We assigned a 
systematic uncertainty associated to this procedure, crji,, of 0.05 
and 0.04 dex for HD 42618 and HD 43587, respectively. 

These various sources of errors (statistical and those arising 
from the mismatch of the reference solar spectrum, the choice 
of the atmosphere models and, in the case of Li, that associ- 
ated to the fitting procedure) can be considered as independent 
and were quadratically summed to obtain the final uncertainties. 
Table |2] shows an example of the error calculations. In the case 
of abundances based on a single line, o-\i„s cannot be determined 
and we conservatively set the total uncertainty to 0.05 dex. 

The results obtained for a given quantity based on the ith line 
list, Xj, were weighted by their total random uncertainties, cr,, to 
obtain the final values, x, provided in the last column of Table [1] 
with the usual formula: 



Uxi/oj) 

m/cr2)' 



(1) 



4. Atmospheric parameters and abundances 

The atmospheric parameters and chemical composition of the 
two Co/?ortargets are provided in Table[T]for the four different 
line lists. The statistical errors are related to the uncertainties in- 
herent to fulfilling the excitation and ionisation equilibrium of 
the iron lines (reff and log g) or constraining the Fe i abundances 
to be independent of the line strength (^). To estimate the uncer- 
tainty in Tetf, for instance, we considered the range over which 
the slope of the relation between the Fe i abundance differences 
and the excitation potential is consistent with zero within the er- 
rors. The effect on the abundances was then examined by altering 
one of the parameters by its uncertainty, while keeping the other 
two fixed. In practice, however, the parameters of the model are 
intimately coupled such that changes in one of these parameters 



The corresponding uncertainty is computed as the quadratic 
sum of the statistical and systematic errors: 



o-( x) 



1 



m/(rj) 



1/2 



(2) 



For HD 42618, we obtain: T^g = 5765+17 K, logg = 
4.48+0.04 dex and [Fe/H] = -0.10+0.02 dex, while HD 43587 
is slightly hotter (in accordance with its spectral type) and metal 
rich: Teff = 5947+17 K, logg = 4.37+0.04 dex and [Fe/H] = 
-0.02+0.02 dex. As shown in FiglH our gravity estimates are 
supported for the two targets by the fact that ionisation balance 
is fulfilled within the errors for three other elements (Si, Ti, and 
Cr). The mean abundance of these elements and iron discussed 
in the following is the weighted (by the inverse variance) av- 
erage of the values for the two ions. The chemical mixture is 
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Table 1. Atmospheric parameters and elemental abundances relative to solar obtained for HD 42618 and HD 43587 using the 
various line lists. 



HD 42618 





Reddy et al. (2003^ 


Chen et al. ll2000ll 


Feltzing &. Gonzalez (|200l!| 


Bensbv et al. (l2003k 


Mean 


Teff [K] 


5745±25 


5765±19 


5765±27 


5775±21 


5765±17 


log^ [cgs] 


4.43+0.07 


4.49±0.06 


4.50±0.06 


4.50±0.05 


4.48±0.04 


f [km s-'] 


1.15±0.06 


0.86±0.06 


1.17±0.06 


0.70±0.05 


0.95±0.04 


[Fe i/H] 


-0.11 ±0.04 (45) 


-0.09±0.03 (82) 


-0.09±0.04 (50) 


-0.11 ±0.03 (120) 


-0.10±0.02 


[Fe ii/H] 


-0.11 ±0.04 (7) 


-0.09±0.04 (7) 


-0.09±0.04 (6) 


-0.11 ±0.04 (25) 


-0.10±0.03 


[Lii/H] 


+0.21±0.07 (1) 


+0.23±0.07 (1) 


+0.23±0.07 (1) 


+0.24±0.07 (1) 


+0.23±0.06° 


[Ci/Fe] 


+0.01±0.03 (2) 








+0.01±0.03 


[Nai/Fe] 


-0.01±0.02 (2) 


-0.02±0.02 (2) 


-0.03 (1) 


+0.00±0.03 (4) 


-0.01 ±0.02 


[Mgi/Fe] 


-0.01 (1) 


+0.04 (1) 




+0.06±0.03 (2) 


+0.05±0.02 


[Ah/Fe] 


+0.00 (1) 


-0.01 (1) 


-0.01 (1) 


+0.01±0.02 (2) 


+0.00±0.02 


[Sii/Fe] 


+0.01±0.03 (6) 


+0.01±0.03 (10) 


+0.00±0.03 (8) 


+0.02±0.03 (16) 


+0.01±0.02 


[Siii/Fe] 


-0.02 (1) 








-0.02±0.05* 


[Si/Fe] 


+0.00 (1) 








+0.00±0.05* 


[Cai/Fe] 


+0.01±0.03 (5) 


+0.01±0.03 (14) 


+0.00±0.04 (7) 


+0.03±0.04 (19) 


+0.01±0.02 


[Scii/Fe] 


+0.00±0.02 (2) 




+0.00±0.02 (2) 




+0.00±0.02 


[Tii/Fe] 


+0.01±0.03 (7) 


+0.02±0.02 (5) 


-0.01±0.05 (7) 


+0.02±0.03 (28) 


+0.02±0.02 


[Tiii/Fe] 








+0.04±0.03(11) 


+0.04±0.03 


[Vi/Fe] 


-0.01±0.03 (6) 


-0.01±0.02 (3) 


+0.00±0.04 (8) 




-0.01 ±0.02 


[Cri/Fe] 


+0.02±0.03 (4) 


+0.00±0.02 (2) 


+0.00±0.02 (5) 


+0.01±0.04 (12) 


+0.00±0.02 


[Crn/Fe] 






+0.03 (1) 


+0.01 ±0.03 (4) 


+0.01±0.02 


[Coi/Fe] 


+0.01±0.04 (3) 




-0.01±0.07 (4) 




+0.00±0.03 


[Nii/Fe] 


-0.02±0.03 (15) 


-0.02±0.02 (16) 


-0.02±0.02 (21) 


-0.01±0.03 (45) 


-0.02±0.02 


[Zni/Fe] 


+0.03±0.03 (2) 






+0.02±0.03 (2) 


+0.03±0.02 


[Sri/Fe] 


+0.00 (1) 








+0.00±0.05* 


[Yii/Fe] 


+0.00±0.03 (4) 








+0.00±0.03 


[Zrii/Fe] 


-0.04(1) 








-0.04±0.05'' 


[Baii/Fe] 


-0.05±0.06 (2) 


+0.00±0.08 (3) 






-0.03±0.05 


[Ceii/Fe] 


+0.03±0.02 (3) 








+0.03±0.02 


[Ndii/Fe] 


+0.03 (1) 








+0.03±0.05* 


[Eu ii/Fe] 


+0.07 (1) 








+0.07±0.05'' 










HD 43587 








Reddy et al. (20031 


Chen et al. (2000) 


Feltzing &l Gonzalez (2001 1 


Bensbv et al. (l2003k 


Mean 


Teff [K] 


5920±22 


5965±21 


5940±27 


5965±26 


5947±17 


log^ [cgs] 


4.29±0.07 


4.42±0.06 


4.39±0.06 


4.38±0.06 


4.37±0.04 


4^[kms-'] 


1.29±0.05 


1.06±0.06 


1.42±0.06 


1.00±0.04 


1.16±0.04 


[Fe i/H] 


-0.02±0.03 (44) 


-0.01 ±0.04 (77) 


-0.02±0.03 (45) 


-0.02±0.04 (118) 


-0.02±0.02 


[Feii/H] 


-0.02±0.04 (6) 


-0.01±0.05 (5) 


-0.02±0.04 (5) 


-0.02±0.05 (23) 


-0.02±0.03 


[Lii/H] 


+ 1.11±0.06(1) 


+ 1.15±0.06(1) 


+ 1.13±0.06(1) 


+ 1.14±0.06 (1) 


+ 1.13±0.05° 


[Ci/Fe] 


+0.00±0.03 (3) 




+0.04 (1) 




+0.01±0.03 


[Nai/Fe] 


+0.00±0.02 (2) 


-0.01±0.02 (2) 


-0.01 (1) 


+0.02±0.04 (4) 


+0.00±0.02 


[Mgi/Fe] 


+0.03±0.06 (2) 


+0.05 (1) 




+0.06±0.02 (2) 


+0.06±0.02 


[Ah/Fe] 


-0.01 (1) 


-0.01 (1) 


-0.01 (1) 


+0.00±0.02 (2) 


+0.00±0.02 


[Sii/Fe] 


+0.00±0.02 (6) 


+0.01±0.03 (10) 


+0.00±0.03 (8) 


+0.02±0.03 (15) 


+0.01±0.02 


[Siii/Fe] 


+0.04 (1) 








+0.04±0.05* 


[Cai/Fe] 


-0.01±0.02 (5) 


+0.02±0.05 (15) 


+0.00±0.04 (9) 


+0.04±0.04 (19) 


+0.00±0.02 


[Scii/Fe] 


-0.02±0.03 (2) 




+0.01±0.03 (2) 




+0.00±0.02 


[Tii/Fe] 


-0.01±0.03 (7) 


+0.00±0.04 (5) 


-0.03±0.05 (8) 


+0.00±0.04 (27) 


-0.01 ±0.02 


[Tiii/Fe] 








+0.02±0.04 (12) 


+0.02±0.04 


[Vi/Fe] 


-0.01±0.04 (6) 


+0.00±0.06 (3) 


+0.00±0.04 (8) 




+0.00±0.03 


[Cri/Fe] 


+0.00±0.03 (4) 


+0.00±0.04 (2) 


+0.00±0.03 (5) 


+0.02±0.03 (12) 


+0.01±0.02 


[Crn/Fe] 






+0.01 (1) 


+0.00±0.02 (4) 


+0.00±0.02 


[Coi/Fe] 


+0.00±0.04 (3) 




-0.04±0.05 (6) 




-0.01 ±0.03 


[Nii/Fe] 


+0.00±0.03 (15) 


+0.00±0.03 (14) 


-0.01 ±0.03 (22) 


+0.00±0.03 (41) 


+0.00±0.02 


[Zni/Fe] 


+0.00±0.02 (2) 






-0.01 ±0.03 (2) 


+0.00±0.02 


[Sri/Fe] 


-0.02(1) 








-0.02±0.05'' 


[Yn/Fe] 


-0.04±0.03 (4) 








-0.04±0.03 


[Zrii/Fe] 


-0.05 (1) 








-0.05±0.05'' 


[Baii/Fe] 


-0.09±0.04 (3) 


-0.07±0.03 (3) 






-0.08±0.03 


[Ceii/Fe] 


-0.04±0.03 (2) 








-0.04±0.03 


[Ndn/Fe] 


-0.02(1) 








-0.02±0.05'' 


[Eu n/Fe] 


-0.11 (1) 








-0.11 ±0.05'' 



Notes. Given are the solar-normalised abundance ratios with respect to Fe (relative to hydrogen for iron and lithium). We use the usual notation: 
[A/B] = [log 6(A) - log e(B)] - [log 6(A) - log 6(B)]© where log 6(A) = 12 + log[A/'(A)/A/'(H)] (N is the number density of the species). Our final 
values for the CoRoT targets are provided in the last column. The number in brackets gives the number of lines used for a given ion. 
° The absolute lithium abundances we obtain are log 6(Li) = +1.15 dex for HD 42618 and +2.05 dex for HD 43587. 
* Fixed error value. 
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Table 2. Computation of the error budget for HD 42618 and the line list of Reddy et al. (I2003I) . 



Statistical eiTors Systematic errors 





O" line 


O" Teff 


0"logg 




O" logg+f/Tcff 


0"Teff+f/logg 


O" Teff+logg/f 


'-'"model 


O" sky Cr (i, 


(X 


Teff [K] 




17 








5 


10 


10 


10 


25 


logg [cgs] 






0.03 




0.04 




0.03 


0.02 


0.01 


0.063 


^[kms-'] 








0.030 


0.025 


0.009 




0.030 


0.006 


0.051 


[Fei/H] 


0.027 


0.015 


0.000 


0.004 


0.011 


0.002 


0.003 


0.002 


0.015 


0.037 


[Feii/H] 


0.014 


0.007 


0.012 


0.006 


0.010 


0.014 


0.002 


0.003 


0.015 


0.031 


[Lii/H] 




0.02 


0.00 


0.00 


0.02 


0.00 


0.01 


0.02 


0.02 0.05 


0.065 


[Ci/Fe] 


0.023 


0.015 


0.005 


0.005 


0.006 


0.003 


0.002 


0.004 


0.003 


0.030 


[Nai/Fe] 


0.007 


0.005 


0.008 


0.003 


0.006 


0.009 


0.000 


0.001 


0.005 


0.018 


[Mgi/Fe] 




0.003 


0.012 


0.000 


0.008 


0.013 


0.005 


0.001 


0.007 


0.050° 


[Ah/Fe] 




0.005 


0.006 


0.005 


0.005 


0.008 


0.002 


0.000 


0.002 


0.050° 


[Sii/Fe] 


0.019 


0.002 


0.005 


0.003 


0.008 


0.006 


0.001 


0.001 


0.002 


0.023 


[Siii/Fe] 




0.023 


0.006 


0.002 


0.009 


0.005 


0.001 


0.003 


0.004 


0.050° 


[Si/Fe] 




0.015 


0.003 


0.005 


0.007 


0.001 


0.001 


0.003 


0.000 


0.050° 


[Cai/Fe] 


0.019 


0.010 


0.009 


0.001 


0.005 


0.010 


0.002 


0.002 


0.003 


0.027 


[Scii/Fe] 


0.011 


0.004 


0.008 


0.004 


0.006 


0.006 


0.009 


0.001 


0.003 


0.020 


[Tii/Fe] 


0.010 


0.015 


0.007 


0.001 


0.002 


0.008 


0.003 


0.003 


0.003 


0.022 


[Vi/Fe] 


0.016 


0.017 


0.006 


0.004 


0.004 


0.008 


0.006 


0.000 


0.002 


0.027 


[Cri/Fe] 


0.014 


0.012 


0.007 


0.001 


0.001 


0.008 


0.001 


0.000 


0.001 


0.022 


[Coi/Fe] 


0.028 


0.007 


0.005 


0.003 


0.002 


0.007 


0.002 


0.001 


0.001 


0.031 


[Nii/Fe] 


0.023 


0.007 


0.005 


0.001 


0.003 


0.006 


0.000 


0.001 


0.001 


0.026 


[Zni/Fe] 


0.018 


0.005 


0.002 


0.001 


0.007 


0.002 


0.004 


0.001 


0.001 


0.021 


[Sri/Fe] 




0.015 


0.008 


0.003 


0.001 


0.008 


0.001 


0.001 


0.004 


0.050° 


[Yii/Fe] 


0.023 


0.003 


0.006 


0.002 


0.005 


0.006 


0.004 


0.002 


0.004 


0.026 


[Zrii/Fe] 




0.004 


0.006 


0.004 


0.005 


0.004 


0.010 


0.002 


0.001 


0.050° 


[Ban/Fe] 


0.055 


0.001 


0.000 


0.010 


0.001 


0.002 


0.008 


0.003 


0.011 


0.058 


[Cen/Fe] 


0.011 


0.001 


0.006 


0.002 


0.008 


0.005 


0.009 


0.004 


0.002 


0.019 


[Ndii/Fe] 




0.001 


0.007 


0.004 


0.010 


0.005 


0.013 


0.003 


0.003 


0.050° 


[Eu ii/Fe] 




0.004 


0.006 


0.005 


0.006 


0.005 


0.010 


0.000 


0.001 


0.050° 



Notes. The line-to-line scatter is given by criinc- Other statistical errors arise from uncertainties in the effective temperature (crxctr), surface gravity 

(o"iogg), or microturbulence (cTf). The changes in, for instance, logg, ^, and abundances associated to varying T^g by its uncertainty are noted as 

o"io6g+f/Tcff • As discussed in Sect|4l the systematic uncertainties are related to the choice of either Kurucz or MARCS model atmospheres (iXniojei), 

the use of a skylight solar spectrum (iXsky), or the subjectivity in estimating the best fit to the lithium feature (cTflt). The last column gives the total 

uncertainty. 

° Fixed error value. 



in both cases close to solar, except perhaps for a slight magne- 
sium excess in the two targets and a slight deficiency of some 
rare, heavy species (e.g., the neutron-capture elements Ba and 
Eu) in HD 43587 (Fig|4|. However, the abundance of these ele- 
ments is based on a few lines that are generally weak and caution 
should be exercised. The star HD 42618 is depleted in lithium to 
similar levels as in the Sun, whereas HD 43587 is dramatically 
more lithium-rich by about one order of magnitude. However, 
its surface Li abundance is by no means particularly high when 
compared to stars with similar temperatures, metallicities, and 
ages (e.g., Sestito & Randich 120051 Ramirez et al. I2012I I. Its 
beryllium abundance is roughly what is expected on the basis of 
its lithium content (Boesgaard et al. l2001l l. The star HD 43587 
is actually part of a quadruple system with not only a pair of 
distant low-mass components (HD 43587B and C; Pravdo et al. 
I2006I I. but also a relatively close (p < 1") and much dimmer M- 
type companion, HD 43587 Ab, following a long {f ~ 30 yr) 
and very eccentric (e ~ 0.8) orbit (Catala et al. 2006] Vogt et al. 
120021 Hart kopf e t al. l20T2l Katoh et al. 2013). As discussed by 
Fuhrmann ( 1201 11 1, neglecting the contamination of the integrated 
spectrum by this cooler component may lead to a small error in 
Teff, logg and [Fe/H] of about 20 K, 0.04 dex and 0.02 dex. We 
therefore cannot rule out that our values are underestimated by 
similar amounts Q 



^ Our spectrum is contaminated by the presence of the companion be- 
cause the fibre entrance aperture projected on the sky of HARPS in the 



5. Reliability of the results 

5. 1 . Fundamental parameters 

As seen in Fig|5] our parameters agree well with the most precise 
results available in the literature. Our Tf^g value for HD 42618 
is identical, within the error bars, with that of Casagrande et 
al. (1201 II ) based on colour indices (5782±54 K). Our value for 
HD 43587 also agrees well with the results of Gray (1994) and 
Bruntt et al. (2004) who used line-depth ratios (5917±21 and 
5923+8 K, respectively). 

A seismic scaling relation using the frequency of maximum 
oscillation power, Vmax, can also be used to precisely infer the 
surface gravity (e.g.. Brown et al. 1991): 



log g = log go + log 



1 , /Teff 

+ Tlog 

^ \ietf,0 



(3) 



This relation is largely insensitive to the effective tempera- 
ture assumed (Areff = 100 K leads to A logg ~ 0.004 dex only 
for Sun-like stars). The seismic and spectroscopic gravities can 
therefore be regarded as being virtually completely independent. 



EGGS mode is 1.4", while interferometric observations - which were 
obtained at the end December 2010 and are therefore nearly contem- 
poraneous with our spectroscopy - give a separation between the two 
components of about 0.88" (Hartkopf et al. i2012t . 
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T" 



T" 



T" 



0.1 - HD 43587 
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Fig. 4. Abundance pattern of HD 42618 (upper panel) and HD 
43587 (bottom panel) relative to the Sun. The abundances of ele- 
ments based on two different ionisation stages (Si i/Si ii, Ti i/Ti ii, 
Cri/Crii, and Fei/Feii) are shown as triangles. For clarity, the 
points have been slightly shifted along the abscissa axis in this 
case. The open symbols denote abundances based on a single 
spectral line, whose uncertainty has been conservatively set to 
0.05 dex. The dashed line shows the mean [Fe/H] abundance. 



Our effective temperatures, Vmax,© = 3100 yuHz, and the pre- 
liminary Vmax values obtained from the analysis of the CoRoT 
data for HD 42618 (3200+100 pHz; Barban et al. 120131) and 
HD 43587 (2300±100 //Hz; Boumier et al., in preparation) were 
assumed to compute the seismic gravities. Propagating the er- 
rors in these quantities leads to very small uncertainties in log g 
(of the order of 0.01-0.02 dex). However, it is very likely that 
these values are unrealistically low and that a systematic error re- 
lated to the scaling relation itself dominates. We conservatively 
adopted a systematic uncertainty of 0.05 dex in the following, 
although the reliability of the relations that link the stellar pa- 
rameters and the seismic observables has been so far validated 
by a number of empirical tests (e.g., Miglio 120121) . In particu- 
lar, there is no evidence for systematic differences from a sta- 
tistical point of view with the gravities obtained using classical 
methods (e.g., the absolute difference for nearby dwarfs with re- 
spect to gravities determined from isochrone fitting is A log g - 
0.006+0.065 dex; Morel & Miglio 2012). We finally obtained 
logg = 4.45+0.06 and 4.32+0.06 dex for HD 42618 and HD 
43587, respectively. As seen in Fig|5] the spectroscopic gravities 
and those based on scaling relations using v^ax agree within the 
errors. A more detailed modelling of the pulsation frequencies 
(additionally using our Tgff as input; see Metcalfe et al. I2012I I 
leads for HD 42618 to a nearly identical seismic logg value 
(Barban et al. l2013l l. The agreement is, however, less gratifying 
for HD 43587, with a spectroscopic estimate higher at the 0.05 
dex level (Boumier et al., in preparation). A similar discrepancy 
is observed when comparing with evolutionary tracks (see be- 
low) and might be related to the presence of the companion. 

A Bayesian estimation through isochrone fitting based on the 
evolutionary tracks of Girardi et al. (I2000 l was also performed 
with PARAM (see da Silva et al. I2006IPI to estimate the funda- 



4.8 



ID 4.6 
-a 



OS 4.4 
o 



4.2 



HD 42618 



T 



CI 1 




5700 5800 5900 

Teff [K] 



6000 



See also |http : //stev . oapd . inaf . it/cgi-bin/param| . 



Fig. 5. Comparison between our Teff and logg values (filled cir- 
cles) and those in the literature (Boesgaard et al. ,2001. Bruntt et 
al. EOO ll da S ilva et al. I MIl Fuhrmann [200 41 EOlI l Gillon & 
Magain'2006'i Gray et al. l2003l Kov tyukh et al. l2004l Mis henina 
et al. 2004 2008, Neves e t al.l2012l Ramirez et al.| 20T3 l Reddy 
et al. 120031 Takeda et al. l2007l Valenti & Fischer 2005 ). The 
values obta ined f or HD 43587 by Bruntt et al. (2004) and Gillon 
& Magain ( I2006I I using different methods are plotted separately. 
Filled triangles: parameters derived from excitation and ionisa- 
tion balance of the iron lines, open triangles: parameters de- 
rived using different techniques. Also shown are the results of 
Casagrande et al. (201 1) for HD 42618 based on colour indices 
(CI 1) and those for HD 43587 of Bruntt et al. (2004, B04) and 
Gray (1994; G04) based on line-depth ratios. The horizontal strip 
corresponds to the 1-cr error bar for the surface gravity deter- 
mined from seismic scaling relations (see Sect l5.U . 



mental parameters. The input values were Teff, [Fe/H], the par- 
allaxes and the magnitudes in the V band (reddening was ne- 
glected). The results are shown in Table |3] The isochrone and 
spectroscopic surface gravities agree within the errors for HD 
42618, while the latter is significantly higher for HD 43587. As 
for the spectroscopic results, the presence of the faint M-type 
companion is also expected to slightly affect these estimates: 
correcting for the light contribution of the secondary (assum- 
ing Ay ~ 4.7 mag; Fuhrmann 120111 1 and revising Teff slightly 
upwai-ds by 20 K and [Fe/H] by 0.02 dex (Sectgll, one ob- 
tains that the star is fainter by about 0.014 mag in the V band, 
slightly younger and more massive (Table|3]l. The small discrep- 
ancy in log g remains when taking into account that our spectro- 
scopic value should also be revised upwards (by about 0.04 dex; 
Sectgll. The radius, R = 1.280+0.032 R0, derived by Thevenin 
et al. (2006) for HD 43587Aa using a calibrated Barnes-Evans 
relation and the previous, slightly lower Hipparcos parallax is 
likely overestimated. On the other hand, its dynamical mass does 
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not provide tight constraints (M - 0.96+0.18 Mq; Catala et al. 
I2OO6I 1. 



Table 3. Fundamental parameters derived using PARAM. 





HD 42618 


HD 43587 


HD 43587 








coiTected 


Age [Gyrs] 


2.17±1.83 


5.69+0.57 


4.97+0.52 


M[Mo] 


0.978±0.026 


1.037+0.016 


1.049+0.016 


[Ro] 


1.00±0.01 


1.18+0.03 


1.15+0.01 


log^ [cgs] 


4.46±0.02 


4.28+0.01 


4.30+0.01 



Notes. The third column gives the PARAM results for HD 43587 Aa 
correcting for the presence of the companion. 



5.2. Elemental abundances 

Our iron abundances are in accordance with those derived by 
previous studies, which range between -0.16 and -0.02 dex for 
HD 42618 and between -0.13 and +0.14 dex for HD 43587 (see 
references in caption of Fig|5]l. On the other hand, empirical cali- 
brations between the strength of some near-infrared features and 
the metallicity yield for HD 43587B a [Fe/H] value fully con- 
sistent with that obtained here for the primary (e.g., Rojas-Ayala 
et al. 1201 0; Terrien et al. l2012l l. Our lithium abundance for HD 
42618 agrees very well with the literature values (Mishenina et 
al. Egos', Ramirez et al.!20T2l Takeda et al. l2007] l. Boesgaard et 
al. ([2001 ) obtained a lithium abundance lower by about 0.2 dex 
in HD 43587 (log e[Li] = 1.84 dex) using parameters very sim- 
ilar to those determined here. A nearly identical value to that of 
Boesgaard et al. ( 1200 11 1 was also obtained by Favata et al. (I1996I ). 
However, their temperature based on photometric data (Tgff - 
5692 K) is significantly underestimated (Fig|5]), and adopting a 
higher value would lead to a Li abundance closer to our estimate. 
Figure |6] shows a comparison between our abundances and the 
results of the most comprehensive studies in the literature for 
these two objects (those of Takeda et al. 120071 were ignored be- 
cause no error bars are provided). Our values suggest an abun- 
dance pattern generally closer to that of the Sun in both cases. 



6. Discussion 

As discussed previously, HD 42618 might be younger than the 
Sun (keeping in mind that ages are very uncertain for such un- 
evolved field stars and that a very young age is not supported 
by the low activity level and lithium abundance), whereas HD 
43587 is perhaps slightly older (Table|3]l. The subsolar metallic- 
ity of HD 42618 considering its relative youth, if confirmed by 
asteroseismology in particular, might have developed because 
it was formed out of more metal-poor interstellar material in 
slightly outer parts of the Galaxy and migrated inwards after- 
wards to its current location in the solar vicinity (e.g., Haywood 
[2008). As a result of the possibly significantly different age of 
HD 42618 compared to the Sun, one might expect subtle dif- 
ferences in its abundance pattern because of the different rela- 
tive proportion of Type la/II supemovae yields and the various 
amounts of material lost by stellar winds from asymptotic giant 
branch (AGB) and massive stars along the history of the Galaxy. 
Tentative (based on data for a dozen stars) relations between the 
stellar ages and the abundances relative to iron have recently 



been proposed for nearby, young dwarfs of the thin disk0 by da 
Silva et al. (I2012I I. However, the abundance pattern of HD 426 1 8 
appears to be very close to solar (among the elements with a de- 
termined uncertainty, magnesium is the most discrepant one but 
by only 2.5 cr; Fig|4|. There is therefore no clear evidence that 
the abundances have been affected by the chemical evolution of 
the interstellar medium. Either the star is not much younger than 
the Sun and/or the differences are too small to be detected. 

Another phenomenon also potentially capable of produc- 
ing slight departures from a scaled-solar chemical mixture is 
the formation of planetary companions. In particular, high- 
accuracy measurements have recently revealed small differences 
as a function of the condensation temperature between the solar 
abundances and those inferred in some solar twins/analogues. 
The different behaviour exhibited by volatile and refractory ele- 
ments is believed to be related to the presence (or lack thereof) of 
terrestrial/gas giant planets (e.g., Melendez et al. 120091 ) . A search 
for planets orbiting our two targets using ultra-precise Doppler 
measurements gave negative results (Howard et al. 120101 Vogt 
et al. I2OO2I 1. As discussed above, the abundance pattern of HD 
42618 is very close to the solar one and no trend between the 
abundance ratios [X/H] and the 50% condensation temperature 
for a solar-photosphere composition gas T^^^^ (Lodders 2003 ) is 
discernible. A similar investigation for HD 43587 reveals that 
none of the slopes between [X/H] and T^^^^^ (either over the 
whole range or over the domain sampled by the refrac- 
tory elements: 7^^^^,,^ ^ 900 K) are statistically significant (see 
FigEJ- Flatter (but still compatible with zero) slopes are obtained 
when the heavy elements (Z > 30) that could bias the results be- 
cause of chemical evolution effects are excluded (see Ramirez et 
al. i201 1) . Flat slopes within the errors have also been reported 
for other solar analogues (e.g., Ramirez et al. 120091 Gonzalez 
Hernandez et al. 120101 Onehag et al. 201 1). The abundances of 
several refractory elements are uncertain in HD 43587. If real, 
however, note that the global depletion relative to solar of these 
elements compared to the volatile ones hinted at by FigQis op- 
posite to the behaviour generally observed in other solar ana- 
logues (e.g., Ramirez et al. l20091 l. 

7. Summary and conclusion 

Our study provides a number of important constraints, which 
will hopefully pave the way for a fruitful modelling of the 
CoRoT data and robust inferences to be drawn about the inter- 
nal structure of these two solar analogues. First, we provided an 
accurate determination of the atmospheric parameters and ele- 
mental abundances through a strictly differential analysis with 
respect to the Sun. While there is a remarkable agreement with 
the surface gravities derived from seismic data or isochrone fit- 
ting for HD 42618, this is less the case for HD 43587. The pres- 

' As described in Morel et al. ( I2003K we used our own HARPS ra- 
dial velocities and the new Hipparcos parallaxes and proper motions 
(van Leeuwen 2007) to compute the space components relative to the 
local standard of rest (LSR). We obtained (f/, V, IV)lsr = (+52.3, -17.9, 
+3.8) and (-29.8, +9.6, -16.0) km s"' for HD 42618 and HD 43587, 
respectively (a right-handed reference system, i.e., with U oriented to- 
wards the Galactic centre was adopted). We then followed Bensby et al. 
( 2004 ) to estimate the probability that the stars belong to a given kine- 
matic component of the Galactic disk. We obtained that these two stars 
are 40 to 50 times more likely to belong to the thin than to the thick 
disk (similar figures have been obtained by Mishenina et al. 120041 and 
Ramirez et al. 120071 1. Adopting another systemic velocity for HD 43587 
Aa+b (the binary system has a semi-amplitude, A" = 4.3 km s"' ; Vogt et 
al. l2002t does not modify our conclusions. 



8 



T. Morel et al.: An abundance study of the solar-analogue Coi?or targets HD 42618 and HD 43587 



0} 



.2 n o a .2 S § 
E-ce>SE.om> 



> 



U 

V 
U 
C 

<a 
c 

XI 



i-, 

a 

"o 

m 

o 



0.2 

0.2 

0.2 


-0.2 

0.2 


-0.2 

0.2 


-0.2 



ECD ^ CD 5 


[Fe/H] 


() 


[C/Fe] 

2S- -J- ■ 




[Na/Fe] 




[Mg/Fe] 

A — 




[Al/Fe] J 
A - 

* : 




[Si/Fe] 




[S/Fe] 
A_ . 




[Ca/Fe] 

zfc-as-i- 




[Sc/Fe] 


-^^^ 


[Ti/Fe] J 




[V/Fe] 




[Cr/Fe] 

■ 




[Co/Fe] 
^$ 




[Ni/Fe] 




[Zn/Fe] J 


L$ 

- ^ 


[Sr/Fe] 

s 




[Y/Fe] 




A 

[Ba/Fe] 




[Ce/Fe] 
25 




[Nd/Fe] J 
1 ^ 

5 i 



Fig. 6. Comparison between our abundances and previous results in the literature (R03: Reddy et al. 120031 VF05: Valenti & Fischer 
M08: Mishenina et al. [20081 GM06: Gillon & Magain l2006t B04: Bruntt et al. l2004l l. The results for HD 42618 and HD 
43587 are shown with circles and triangles, respectively. Our abundances based on a single spectral line have their uncertainty 
conservatively set to 0.05 dex. 



0.05 


-0.05 

-0.1 

-0.15 

-0.2 
0.05 



-0.05 
-0.1 
-0.15 



HD 42618 



slope = -0.34±1.52 

slope = +1.57±2,97 



S Zn Na 



H \ \ 1 — 

HD 43587 



H 1 h 



- - slope = -1.32±1.58 
■■ slope = -3,13±2.98 




500 



1000 



1500 



Fig. 7. Abundances as a function of the condensation tempera- 
ture. Same symbols as in Fig|4] The fits weighted by the inverse 
variance, either taking into account all the elements or only those 
with r^" , > 900 K, are shown as dashed and dotted lines, re- 

cond ' ' 

spectively. The slopes (in units of 10"^ dex K ') are indicated. 



ence of a nearby, low-mass companion may contribute to these 
discrepancies, which remain small and comparable to the uncer- 
tainties, however. The abundance pattern of most metals rela- 
tive to iron appears to be very close to solar for most chemical 
species (unlike what has been found for other solar-like CoRoT 
targets; Bruntt '2009), except for magnesium and a few trace el- 
ements in HD 43587, whose abundances are generally uncer- 
tain. On the other hand, solar abundances within the errors have 
been reported in these two targets for four elements not anal- 



ysed here, namely O, K, Mn, and Cu (Bru ntt et al. 120041 Gillon 
& Magain '2006'; Mishenina et al. l2U08l Ramirez et al. l20T3t 
Reddy et al. 2003). Assuming that the extent of microscopic 
diffusion is similar to that in the Sun, it is thus reasonable to 
assume a scaled-solar composition when modelling the CoRoT 
data. Although we are unable to precisely determine their rota- 
tion rate, these stars are unlikely to be fast rotators unless they 
are seen from an unfavourable inclination angle (i.e., nearly pole 
on). Although not decisive in itself, their low activity levels seem 
to argue against this possibility. Generally speaking, the level of 
lithium depletion in the interior of main-sequence stars tends to 
increase with stellar metallicity and age, but to decrease at higher 
masses (e.g., Li et al. 2012 ). Detailed modelling is required to in- 
vestigate whether the different Li surface abundances observed 
in our targets and in the Sun simply arise from their different 
fundamental parameters or instead reveal internal conditions di- 
versely favourable to the destruction of this fragile element. A 
determination of the rotation properties and depth of the con- 
vective envelope from asteroseismology would be particularly 
valuable in this respect. 
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